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Abstract 



< 

In this article we examine the densities of a product and a ratio of two real positive definite matrix- 

variate random variables X\ and X2, which arc statistically independently distributed, and we consider 

i i 1 1 

the density of the product U\ = X 2 2 X\X 2 as well as the density of the ratio U2 = X 2 Xi X. 2 2 . We define 

matrix- variate Kober fractional integral operators of the first and second kinds from a statistical perspective, 

making use of the derivation in the predecessor of this paper for the scalar variable case, by deriving the 

densities of products and ratios where one variable has a matrix- variate type-1 beta density and the other 

variable has an arbitrary density. Various types of generalizations are considered, by using pathway models, 

by appending matrix variate hypergeometric series etc. During this process matrix-variate Saigo operator 

and other operators are also defined and properties studied. 

o\ 
cn 

All the matrices appearing in this article are p x p real positive definite. Corresponding results for 
hermitian cases can be obtained but not given in this article. The following standard notations will be used. 
A > O means the p x p real matrix is symmetric, X — X' , and further, it is positive definite. \A\ means the 
determinant of A, tr(A) = trace of A, dX will stand for the wedge product of differentials in any matrix X. 
If X is p x q with X = (xij ) then 



H 

■ dX = dxxx A dxi2 A ... A dx pq for a general matrix (1.1) 

= Y[ Adxy = J^J Adxij when X = X' or when X is symmetric. (1-2) 

i>j j>i 

Also J x f(X)dX will mean the integrals over all X (need not be symmetric or even square) of a real-valued 
scalar function f(X) of matrix argument X . In the same format 

f(Y)dY= [ f(Y)dY 

i JO<A<Y<B 

will mean the integration of a real-valued scalar function of the real positive definite p x p matrix Y over 
the space of positive definite matrices such that A>0,X>0,X — A>0,B — X>0 where A and B are 
p x p constant matrices. The notation will then imply that if O < X < I then all eigenvalues of X are in 
the open interval (0, 1). We will need some Jacobians of matrix transformations in this paper. For further 
results on Jacobians, see Mathai (1997). 
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Y = AX A', X = X', | A\ ^ => dY = \A\ p+1 dX (1.3) 

! , / |X|- 2 PdX for a general X . , 

Y-X -{ ]xr ( P +t) dx iorX = x >. ( 1A > 

We will denote the unique positive definite square root of a positive definite matrix A by A*. The following 
standard property will be used very often in this article. For px p nonsingular matrices A and B 

\I ± AB\ = \I± BA\ = \A\ {A' 1 ±B\ = \B\ \B~ l ± A\ (when nonsingular) 

\I±AB\ = \I±A*BA*\ = \I±B$AB*\ (when positive definite). (1.5) 



The real matrix- variate gamma function, denoted by Y p {a) is defined as follows which has an integral 

3—] 

T" 



representation when 5R(a) > ^— - . 



T p (a) = n^Tiana ±)...r(a ^), H(a) > ^ 

= / |X| Q -4V tr ( x )dX, R(a) > ^— -. (1.6) 
Jx>o 2 

The real matrix- variate type-1 beta density for the pxp positive definite matrix X\, with parameters a and 
/3 and denoted by /i(Xi), is defined as follows: 

^ = r^rl ^-^ 1 -*^' 0< Xl <i 

T p {a)T p {fj) 

for 5R(a) > SR(/3) > ^ and /(Xi) = 0, / 2 (Fi) = elsewhere. Type-1 and Type-2 beta integrals and 
beta functions are defined and denoted as follows for 5i(a) > 3?(/3) > 

B p (a,0) = Tp{a)TpiP) 



-L 



X\ a -^\I - X\"-^dX (type-1) 



o<x<i 



= [ lYf-^lI-Y^-^dY (type-1) 
Jo<y<i 

= / l^r-^IZ + CTl-^+^d^ (type-2) 
Ju>o 

= [ {Vf-^ \I + V\- {a+f}) dV (type-2) (1.7) 
2. Kober Operator of the Second Kind for the Real Matrix-variate Case 

Definition 2.1. We will define and denote Kober operator of the second kind for the real matrix-variate 
case as follows: 

K<x a f(X) = / \T-Xr P -¥\T\-<- a f(T)dT,Vt(a) > ^ (2-1) 

i JT>X z 



Consider two il matrix- variate random variables X\ and X2, independently distributed, where X\ 

0+1 
2 



has a real matrix-variate type-1 beta density fi(Xi) with parameters (( + a), that is, 



h{Xl) = r P (C + <* + p + 1 ) c _ ^ 

for 3?(C) > — l,3?(a) > and /i(Xi) = elsewhere. Let X2 have an arbitrary density /(Xj). Then the 

11. 

joint density of X\ and X2 is f\{X{)f{X2). Let us consider the transformation U = X%X\X% , X 2 = so 
that X 2 = V,X 1 = V-^UV-? and the Jacobian is given by cLYi A dX 2 = {V^^dU A dV. If the joint 
density is denoted by /([/, V) then 

m v)du a dv = lf^"l%\ \v-iuv-i\< 

x IJ-V-il/V-'l^^/^lVI-^dt/AdV: 

Therefore the marginal density of U, denoted by g(U), is available by integrating out V from f(U, V). That 
is, 



g(u) 



[ h{v-*uv-*)f{v)\v\-*&v 

JV 

r P (C + a + H±i) f i 



L p (C+^i) iy>[,r p (a) 



x |^|«|vr c |vr a+£ * 1 |y-Cr| a - ,! i ± /(V)dV 

r p (a + c + y) 
r P (C + 2 } 1 ) " 

Hence we have the following theorem: 

Theorem 2.1. When X\ and X2 are independently distributed p x p positive definite real matrix random 

II _i _i 

variables and when X 2 = V and U = X^X X X^ or X X = V 2 JJV 2 and when X\ has a real matrix-variate 

type-1 beta distribution with the parameters (( + ^^,a) and if g(U) denotes the density of U then 



r P (C + ^) 
r p (« + c + ^) ; 



-a(U) = K^f{U) (2.2) 



2 

is Kober fractional integral operator of the second kind for the real matrix-variate case. 

As a special case of (2.2), or independently, we can derive a result for the right sided Weyl operator for the 
real matrix-variate case. Let the right sided Weyl fractional integral operator be denoted by (xW^ a f){X) 
with infinity here signifying that T — X is positive definite. 

Theorem 2.2. Let Xi,X2,U,V be as defined in Theorem 2.1. Let X\ have a type-1 beta density with the 
parameters (^-,a). Let the arbitrary density of X 2 be denoted by / 2 (X 2 ) = |A" 2 | Q /(A" 2 ), where f(X 2 ) is 
arbitrary. Let the density of U be again denoted by g{U). Then 



( x w- a f)(x) = -L / \v- up-*? f{v)dv 

r (a) Jv>u 



^T-g(U), R(a) > (2.3) 
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2.1. A pathway generalization of Kober operator of the second kind in the matrix case 



A pathway generalization, parallel to the results in the scalar case can be considered. In the pathway 
case when generalization to matrices is considered we take 5=1. For a general 5, there will be problems 
with Jacobians of transformations for X even if X > O and 5 > 0, see for example Mathai (1997) for the 
case 5 = 2 and when X = X' to see the type of complications. Hence we take the case 5=1 only. Let X\ 
have a pathway density 

h(X 1 ) = C 1 \X 1 y\I-a(l-q)X\^ (2.4) 



for I — o(l — q)X >O,q<l,r]>0,a>0 where C\ can be seen to be the following: 

_ _ Ki-g)F 7+£ ^r p ( 7 + T iL + (p + i)) 
' r p(7 + ^)r p(T ^ + H±i) 



Let X 2 have an arbitrary density f(X 2 ) and let X\ and X 2 be statistically independently distributed. Let 

U = XlX x Xl,X 2 = V or Xi = V'sUV-i and the Jacobian is l^l -2 ^. Let g(U) be the density of U. 
Then, going through the earlier steps we have the following: 

g(U) = Ci\U\* f [ {Vr'-^+^lV -a(l-q)U\^f(V)dV. (2.6) 

JV>a(l-q)U 



Then 



r P (7 + — meo = 5 p+1 i^r 

i PVl- 9 "+" 2 / 

/ |T/-a(l- g )[/|T^|y|-T-(T^+ £ r i )/(]/)dT/ 



y>a(l-g)?7 

= K u5 +E ^ f(U) (2-7) 



where Ky^~* 2 f(U) can be called the generalized pathway Kober operator of the second kind in the real 
matrix- variate case. When the pathway parameter q varies from — oo to 1 it provides a pathway or a class 
of operators and all these operators in this pathway class will eventually go to the exponential form. For 
a = 1, q = 0, j- 2 ^ = a — 2±I and 7 = ( we have 

K 1 u ^ +V ^ L f{U) = K^ a f{U) (2.8) 

the Kober operator of the second kind as a constant multiple of the density of the product of two matrix- 
variate independently distributed random variables. Note that when q — > 1_ we can evaluate the limit of 
g{U) by using the following lemmas: 

Lemma 2.1. 



= («£ 

r p ( 7 + 2±i) u 

Proof: Open up each T p (-) in C\ of (2.5) in terms of ordinary gamma functions. Then use the following 
asymptotic approximation for gamma functions. For \z\ — > 00 and 7 a bounded quantity 

r(z + 7 ) « T^z^-'-Se- 2 . (m) 

This is the first term in the asymptotic series. This term is also known as Stirling's approximation. When 
q — > 1_ we have — > 00 and hence take \z\ as and expand all gammas by using Stirling's approximation 
to see that C\ reduces to (i) above. 
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Lemma 2.2. 

lim 17 - o(l - q)X\& = c- a " tl(x K (Hi) 
i -hi- 
proof: Writing the determinant in terms of eigenvalues we have 

\I - a(l - q)X\^ = JJ(1 - a(l - q)\j)^ (iv) 

3 = 1 

where Ai, A p are the eigenvalues of X. Now 

lim (1 — a(l — q)\j) T ^ = e~ ari ^ j . (v) 

Hcncc 

p 

lim |/ - o(l - q)X\^ = Y[ c~ Wn Aj = e _< "> tr W 
9 ~* 1 ~ 3=1 

which establishes (hi). Now by using Lemmas 2.1 and 2.2 we have 

lim g(U) = {a7ir+P ~T \UV 

x / \ V \-y-^e- a, i^ v ~ iuv ~^dV. (2.9) 
Jv>o 

This is the limiting form of the pathway Kober operator of the second kind in this class of pathway operators 
of the second kind. 

In the pathway generalization, one can also replace the parameter a with a constant positive definite 
matrix A. In this case the model will be written as 

f\{Xi) = dWlXiWl - (1 - q)A*XiA*\&* (2.10) 
for q < 1, A > O, X\ > O, I — (1 — q)A^X\A^ > O. The pathway parameter is still q. In this case 



r(M (i- g r +i ^l^ + ^r p(7 + 7 ^ + ( P + i)) 

CM) = f^Tfj^ • (2 - n) 

Then g(U) of (2.6) goes to the following form, denoted by gA(U) 

g A (U) = C 1 (A)\A\^\U\^ 

x f \V^A- X V^ - (l-g)?7|T^|y|-^-(T^+4 i )/(y)dT/ (2.12) 
Jv*>o 

where 

V* =viA~ 1 vi - (l-q)U. 
Then one can define a pathway generalized Kober operator of the second kind as 

Kut^fV) = r p(7 + P ~^)9a{U) 



_ (i - ,r+^ i Ar^ + ^ Tph + + (p + 1)) 

T { _n i p±1\ 

X / |^3A- 1 y^-(l- (? )C/|T^|l/|-T-(T^ + £ T i )/(-l/)dT/ (2.13) 

Jv*>o 



In this case, as q — > 1_ we have 
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lim g A (U) = i-i \UV 

9-1- r p ( 7 +H+i) 1 1 



x 



f \ V \-l- P -¥ e -ntr(Ahv-huV-hAh) f{ y )d ^ (2W) 
JV>A 



3. M-transforms of Kober Operator of the Second Kind 

The generalized matrix transform or M-transform is defined and illustrated in Mathai (1997). The 
M-transform of Kober operator of the second kind is the following: 

Theorem 3.1. For the Kober operator of the second kind defined in (2.1) the M-transform with parameter 
s is given by 

M{K^ a f(X);s}= [ \X\'-^[f ^\T-X\ a - t ¥\T\-<- a f{T)dT\dX 
Jx>o Jt>x 1 p\ a ) 

-WTJTTf^ *(C + -)>^,*(«)>^ (3-D 

where f*(s) is the M-transform of f(X). 
Proof: Interchanging the integral we have 

M{K^ a f(X);s}= [ \T\-<- a f(T){-^ [ |^|C+— ^ | T _ X \ a -^dX]dT. 

JT>0 L P\ a ) JO<X<T 

Note that 

\T-X\ = \T\ \I -T~iXT-i\, Y = T~^XT~i dY = \T\-^dX. 



\X\t- E ¥\T-X\ a - E ¥dX=\T\^+ s - E ¥ / \Y\<+ s - P -¥\I-Y\>*- P -¥dY 



Hence 

■ \ / Ji 
We can evaluate the Y~-integral by using real matrix- variate type-1 beta integral. 



/ 

Jo 



\Y\ c+s -**-\I - Y\ a -' L s-dY 



r p (( + S )T p (a) 



lo<Y<i r p (a + c + s) 

for K(a) > ^i,5R(C + s) > 2^-. Now the T-integral gives 

f \T\ s - E ¥f(T)dT=f*(s) 
Jt>o 

where f*(s) is the M-transform of f(X). Hence (3.1) follows. Note that for p = 1 the result agrees with 
that in the scalar case, which is available in the literature, see for example Mathai and Haubold (2008). 

From (3.1) for ( = and 31(a) > we have the special case of the Kober operator of the second kind 

K°x a f{X) = ^r-J \T-Xr^\T\- a f(T)dT. (3.2) 
1 \ a ) Jt>x 

But the right side of (3.2) is Weyl fractional integral of order a in the matrix case, xW^ ) a f(X), except for 
the factor |T|~ Q . The Weyl integral in the matrix case is 

xW~ a f(X) = ^r-J \T-Xr^f(T)dT^(a) > (3.3) 
1 p\ a ) Jt>x z 



Hence we have the following corollary. 

Corollary 3.1.1. The M-transform of the right sided Weyl operator in the real matrix case is given by 

M{ x w- a \x\-«f(x) ]S } = rffi^ rOO (3-4) 

for sft(s) > 2zi, K(a) > ^ where f*(s) is the M-transform of f(X). 

The proof is parallel to that in Theorem 3.1. Let us see whether a Mellin convolution type formula holds for 
Kober operator of the second kind in the matrix case. Let 

g(U)= [ \V\~ P -^ h{V-^UV-^)f 2 {V)dV (3.5) 
Jv 

where fi(Xi) is a type-1 matrix-variate beta density with parameters (C + , a). That is, 



L p (a)L p (C+2±i 

for K(a) > ^-,^(0 > -1 and /i(Xi) = elsewhere. Substituting (3.6) in (3.5) we have 



Lp(a)L p (C+ ^i-) 



x iJ-y-^c/y-sp-^/^dy 



r p (a) 

= K^ a f(U) (3.7) 
which is the Kober operator of the second kind. Hence we have the following theorem: 

Theorem 3.2. Kober operator of the second kind with real matrix argument can also be represented as a 
Mellin convolution type formula 



K C x a f(X) = f \V\-^ fi{V-*XV-t)f2(V)dV 
Jv 



where f\(X\) is a type-1 beta density with parameters (C + ^-,a) and /bOO is an arbitrary function or 
arbitrary density if the Kober operator is to be taken as a constant multiple of a statistical density. 

4. Generalization in Terms of Hypergeometric Series for Kober Operator of the Second Kind in the Real 
Matrix Case 

For introducing hypergeometric series of matrix argument we will need the definitions, notation and lem- 
mas. Hypergeometric functions of matrix argument are defined in terms of matrix-variate Laplace transforms, 
M-transforms and zonal polynomials. Explicit series form for all cases is available through the definition in 
terms of zonal polynomials and hence we will define in terms of zonal polynomials. 

r F s (Z) = r F s (ai, a r ; bi, b s ; Z) 

\ ^ \ ^ (a)K — {ar)K C K {Z) . . 

where K = (fci, k p ), k\ + ... + k p = k is a partition of k = 0, 1, 2, ... 
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(a) K = f[(a- ] -^-)k 3 , (b) kj = b(b + 1)...(6 + kj - 1), (b) = 1,6^0 (4.2) 
j=i 

and Ck{Z) is a zonal polynomial of order k and Z is a p x p matrix. The series is defined for the real 
and complex matrices. Zonal polynomials are certain symmetric functions of the eigenvalues of Z. In our 
discussions, Z will be real and positive definite. For more details about zonal polynomials see Mathai, Provost 
and Hayakawa (1995). The following basic results are needed in our discussions. A standard notation in this 
area is 

r p (a,K)=T p (a)(a) K . (4.3) 
The following basic results are needed in our discussion. 

Lemma 4.1. 



Xr^ll Xf- P -^C K {TX)AX = 1 ~ ,L a C K (T) (4.4) 



T p {a,K)T p {p) i 



T p (a + p,K) 
for^{a) > Efi.ROS) > 

Lemma 4.2. For R(a) > *fi, A > O, S > O 

I \Sr^C K (ZS)dS = T : { " ,K J p f~tl \ArC K (ZA). (4.5) 
Jo<s<a T p (a+ ! ^-,K) 

Let us assume that all the parameters ai, a r , &i, b s are real and positive and let the argument 
matrices be p x p and positive definite. For A > O, let the density of X\ be 

= - r F s ( ai , a r ; h, b s ; AX 1 )\X 1 \<\I - X^- 2 ^ 

- - E E - ^r-** 1 (4.6) 

where the normalizing constant c/ is available by integrating out term by term with the help of Lemma 4.1. 
It will be available in terms of a r+iFj+i. Let /(X2) be an arbitrary density. As before, let Ai = ^"^{TV - ^, 
then denoting the density of U again by g(U) we have 



9{U) 



f .h{V-^UV-^).f{V)\V\^dV 

J V 

1 T p (a + C+^-) ^y^y-v (ai)j<:-(a r )if 1 



c /r p (C + ^)r p (a) v ^ ^ (6i)k-(6 s )k fc! 

x / |V , -'^V'-'| c |/-V'-'^V , -'| a - £ * 1 |V'|- £ * 1 Cjf(ilV , -'I7V'-')/(V')dV' (4.7) 

This is the generalization of a constant times the Kober operator of the second kind in the matrix case. For 
r F s = 2-F1 one has the matrix-variate generalization of a constant times the Saigo operator of the second 
kind in the real matrix-variate case. 

5. Kober Fractional Integral Operators of the First Kind in the Matrix Case 

Definition 5.1. Wc will give the following definition and notation for Kober fractional integral operator 
of the first kind in the real matrix-variate case: 

T (,a »/^n |A| a f T r,™_£±i, 



/r/< V) = ffi-^ / \X-Vr^\Vff(V)dV (5.1) 
ip(aj Jv<x 
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for sft(C) > -l,K(a) > 

Theorem 5.1. For $l(a) > > —1 ifte M-transform, with parameter s, of Kober operator of the 

first kind in the real matrix-variate case, is given by 



M{4 a f(X);s}= [ \X\'-^[^-^- [ \X-V\ a - E ^\V\if(V)dV}dX 
Jx>o L p{ a ) Jv<x 

r (t + £±± - a) n _ 1 

= 2 > p, s R / s < R ,£ + x) R / ) > 

r P (« + C + ^ - *) 2 

w/iere /*(s) «s t/ie M-transform of f(X). 

Proof: Integrating out X first we have the X-integral 

f | X | S -C- a -£+i| X _ y |a-£+l dx= /* | y + y | S - C - a -£+i | y | a _£+i dy y = ^ _ y 

Jx>v Jy>o 



>Y>0 

Put Z = r^yy-i => dZ = \V\-^-dY. Then the X-integral is 



= | y | 5 - C -a-£±I f | /+F _i yF _i| 5 _ C _ a _£+i| y | a _£+i dy 

Jy>o 



|Zr" £ i i |7 + Z|-( E ? i +^- s )dZ 

z ( ) 

,_ c _i3i r p (a)r p (g±i + c-«) 
r p (£±i + a + c- s ) 



for 3?(a) > ^2^", 3?(C — s) > — 1 by evaluating the integral by using a type-2 matrix-variate beta integral in 
the real case. Now, the F-integral becomes 

/ \vr p -^f(v)dv = r( S ). 

Jv>o 



r p (^ + a + c-sy 



M{4 a f(X); s} = Z/^, "» - f*(s) (5.3) 



Hence 

il/f/K' ^ c\ - I 

(P±l 

for 5R(a) > ^^(C - s) > -1. Note that for C = 0, 

I° x a f(X) = \X\-%D x a f(X) (5.4) 

where oF> x a is the left sided Riemann-Liouville fractional integral for the matrix-variate case. Note that for 
the scalar case, for p = 1, 

^■ a /(»)^> = r ( i + a + c- fl ) (5 ' 5) 

for 3?(a) > 0, 3ft (£ — s) > — 1 agreeing with the corresponding Mcllin transform in the scalar case. 
Corollary 5.1.1. The M-transform of \X\- a D x a f (X) is given by 



for K(a) > ^.^(s) < 1. 



The proof is parallel to that in Theorem 5.1. 

Let us treat a Kober operator of the first kind as a statistical density. Let X 2 have an arbitrary real 
matrix- variate density f(X 2 ) and X\ have a real matrix- variate type-1 beta density with parameters ((, a). 
That is, 

= fW^ if-^l* W^'OK X 1 < I (5.7) 

r p(C)r P (a) 

for Sft(£) > £ rp, ^(a) > £ rp and fi(Xi) = elsewhere. Let X\ and X 2 be statistically independently 
distributed. Consider the transformation X 2 = V,X\ = V^U~ 1 V% . The Jacobian is given by 

dX 1 AdX 2 = iV^pl-^+^dU AdV. 
The marginal density of U, denoted by g r (U) where r designates that it is coming from a ratio, is given by 



r p(C)r P (a) Jv 

x\I- Viu^vi^-^ 1 fiV^Vl^pl-^+^dV 

= r P ^r rV '" C "° / \U-Vr P -^\V\<f(V)dV. 
T p (QT p (a) J v<u 



Therefore 



rp(0 :9r(U) = l -^Ll / \V\<f(V)dV 



r p (( + ay' y ' T p {a) 



v<u 



= I C x a .f(X). (5.8) 

This is Kober operator of the first kind in the real matrix- variate case and it can be considered as a constant 
multiple of a real matrix- variate statistical density. 

One can also consider a pathway extension for the real matrix-variate Kober operator of the first kind. 
5.2. Pathway Extension of Kober Operator of the First Kind in the Matrix Case 

Consider the following pathway modified form of the density for X\. That is, 

h(X 1 )=C 2 \X 1 p- E i 1 \I-a(l-q)X 1 \^,I-a(l-q)X 1 > O (5.9) 
for q < 1, a > 0, r\ > where 

C 2 = — — P ' . 1-9 2 —. (5.10) 

r p(T^ + l T i ) r p(7) 

Consider the same type of transformation as before: X 2 = V,X\ = vi[J~ 1 Vi . The marginal density of U, 
denoted by g p (U), is given by 



o{U) = C 2 [ \V?U- 1 V?\~<- E i 1 \I-a(l-q)V?U- 1 V?\ 
Jv 



9 P { 

iv 

f{V)\V\^\U\- (p+1) dV (5.11) 



= C 2 |C/r 7 - ( T^+ £ * i ) f \U -a(l-q)V\^\Vpf(V)dV. 

JU>a(l-q)V 

Then 

[a(l - <z)]^L p ( 7 + -IL- + H±i) 
r P (l)9 P (U) = ^ ; T T p+T 2 —\U\-^r- q ^) 



[ \U -a(l-q)V\^\V\y(V)dV (5.12) 

JU>a(l-q)V 



l(l-q) 
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The right side of (5.12) is the pathway extension of Kobcr operator of the first kind. The right side divided 
by r p (7) is also a statistical density of a type of ratio of independently distributed matrix- variate random 
variables. 

Note that for a = 1, q = 0, ^ + 2±1 = a, (5.12) reduces to the special case (5.7) for 7 = (. Thus, (5.12) 
describes a vast family of operators which can all be considered as generalizations of the Kober operator of 
the first kind in the real matrix- variate case. The limiting form when q — > 1_ is available from the structure 
in (5.11). Note that 



lim \I-a(l- q)V^XJ- x V^ 



I , 3 

1-9 



-ar] tr(V3E/ _1 V5) 



q-H 



(5.13) 



Hence 



lim g p (U) = ( lim C 2 ) f \U\ 
q^ 1 - Jv>o 

ar, tr(yi[/-Vi) f (^V)dV 



P+1 

T-T-|y|7 



x e 



(5.14) 



where 



That is, 



r p ( 7 ) 



9^1- r p (7) 



|c/|- 



/ \V\ 7 e' an tl(viu lyi) f(V)dV. 
Jv>o 



(5.15) 



In this case also one can replace the parameter a in /i(Ai) by a constant positive definite matrix A > O. 
Then fi(Xi) denoted by /ia(Ai) can be written as 



where 



f 1A (Xi) = C 2 (A)\X^- — \I-(l-q)A^X 1 A^\- 

(l-g)P7|^|7r p ( 7+T a_ + E±l) 



C 2 {A) 



r P (7)r P (^- + ^) 



Then the density of U = V? U 1 V? , denoted by gA(U), is given by 



<u(£0=C 2 (i4)|tf| 



-7-( T ^+ £ f i ) 



</ L |Fn^-(l-g)ViA^i|^«/(V)dV: 

Jc/>(i- 9 )y 2 2 



(5.16) 
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